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ABSTRACT 



A study of the parameters involved in electron 
scattering experiments with gaseous targets at the Naval Post- 
graduate School Linear Accelerator was conducted. The para- 
meters included gas densities, effective target thicknesses, 
and radiative corrections to experimental cross-section cal- 
culations. To compare the NPGLINAC with other accelerators, 
an electron scattering experiment was performed, using a 
hydrogen-helium gas mixture target, at values of from 
0.10 to 0.55 F~2. The resulting charge form factor for he- 
lixm was fairly consistent with previous experimental models; 
however, it was concluded that experiments using separate 
gas targets yield greater accuracy than those using a gas 
mixture . 
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I . INTRODUCTION 



A. EXPERIMENTAL OBJECTIVES 

The electron linear accelerator at the Naval Postgraduate 
School (NPGLINAC) has been in use since 1966, and has been 
described in theses by Barnett and Cunneen [Ref.lJ and 
Midgarden [Ref. 2]. Energy loss and nuclear structure exper- 
iments have been conducted using solid targets. In 1970, the 
scope of experimentation was broadened to include gas targets 
in addition to solid targets. While the technique of elec- 
tron scattering measurements remained the same, there were 
several new problems to be considered v;hen using gas targets 
instead of solid targets. Some of these problems included 
the following mechanical considerations : 

(a) Containers had to be designed to hold gases under 
high pressures. 

(b) A method of cooling the target area and keeping it 
at a constant temperature had to be designed and 
installed. 

(c) A method of feeding and evacuating the gas target 
chambers had to be devised, along with a method to 
measure the pressure inside the chamber. 

The problem associated with high pressures and low tempera- 
tures arose because it was necessary to maintain the gases 
at high densities in order to obtain reasonable counting 
rates during experiments. Gases at atmospheric pressure and 
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temperature produced counting rates that were very low and 
much more time would have been required for a run. 

These mechanical problems were solved by Professors 
Bumiller and Buskirk, Captain Louis Gaby, and the NPGLINAC 
technicians. A diagram depicting one of the gas target 
chamber arrangements is given in Figure 1. 

Apart from the mechanical problems encountered when gas 
targets were introduced, several experimental problems also 
were evident. The experimental considerations included the 
following ; 

(a) Effects of the target chamber windows on electron 
scattering experiments had to be determined. 

(b) The effective target thicknesses of gases had to be 
determined; i.e., only electrons scattered from a 
small part of the total volume of gas in the target 
chamber made their way into the spectrometer 
entrance . 

(c) The densities of the gases at different pressures 
had to be determined for use in cross section 
calculations . 

These and other experimental considerations have been applied 
to the gas target experiments conducted at the Naval Post- 
graduate School. So far pure deuterium and hydrogen-deuter- 
ium mixtures have been studied by Gaby [Ref. 3] and Mader 
[Ref. 4]. Hydrogen-helium mixtures have also been studied and 
are reported in this thesis. 
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Figure 1 . Gas Target. 
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The purpose of this paper is to provide guidelines for 
use in further studies of gas targets at NPGLINAC. Data 
from deuterium and helium experiments were used in preparation 
of this thesis. 

B. EXPERIMENTAL CALCULATIONS 

One of the most important considerations in a cross- 
section experiment such as a gas target study is the evalua- 
tion of radiative corrections. They are of three types: the 
ionization correction, the Bethe-Heitler correction, and the 
Schwinger correction. These corrections are described in 
detail by Gordon [Ref. 5]. The experimental problems of tar- 
get chamber window thicknesses, effective gas target thick- 
nesses, and gas densities all appear in the Bethe-Heitler 
correction. It was therefore considered constructive to study 
these parameters and report their significance on experimen- 
tal outcomes. Likewise, the Schwinger correction was studied 
for similar reasons. The Schwinger correction has many 
forms and it was necessary to decide which form was to be 
used by students. Conformity is important because typical 
Schwinger corrections are of the order of 10-20%. 

In addition to the study of the effects of experimental 
parameters on radiative corrections, a cross-section exper- 
iment was performed using a hydrogen-helium mixture as a 
target. From the relative cross sections, a calculation of 
the helium charge form factor as a function of the 
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four-momentvim transfer squared v;as performed. The results 
were compared to those obtained at Stanford in 1967 lRef.6] 
and to the Gaussian model of the form factor. 

II. RADIATIVE CORRECTIONS 

A. BETHE-HEITLER CORRECTION 

The Bethe-Heitler radiative correction was developed to 
account for part of the tail below a peak in an electron 
scattering experiment. After incident electrons have scat- 
tered off a target nucleus, they lose energy in the form of 
bremsstrahlung as they pass other nuclei and the atomic elec- 
trons of the target nucleus itself. The less energetic 
electrons are counted at much lower energies than the peak 
energy, theoretically, all the way to zero energy. It is 
not possible to set the spectrometer all the way to zero 
energy because of the length of time involved (the ten- 
channel system can only monitor an energy range of about 
three Mev per setting) and the background radiation levels. 
The Bethe-Heitler correction is an estimate of the n\amber 
of electrons that should be in the radiative tail due to 
bremsstrahlung . 

Of the different forms of the Bethe-Heitler correction 
term, the one used was developed by Tsai [Ref. 7] and appears 
in a thesis by Gordon [Ref. 5]. 
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( 2 - 1 ) 



= I K Tiw + 1/2 bt] In (E^/n^AE) + 

Tf^ + 1/2 bt] In (E 3 /AE)] 

where the b's are functions of the atomic number g only, 

Tj_^ = entrance windov/ thickness 
Tfv; = final window thickness 
t = effective target thickness 
E^ = incident electron energy 
E3 = elastic peak energy 

AE = E 3 - (lower limit of electron spectrum) 

and ri =• recoil factor = 1 + 2E-^ SIN^ 9 

Mc2 2 

where M = mass of target nucleus 
and 6 = scattering angle. 

The Bethe-Heitler correction is applied to the differen- 
tial cross section in the form Kg = e~'^B. in a ratio experi- 
ment using a two-gas mixture, a relative cross section is 
determined, and the ratio of Kg for the two gases is enter- 
ed. The ratio is very nearly equal to one in most cases. 

The Bethe-Heitler term becomes much more significant in a 
pure gas experiment. For the hydrogen-helium runs the ratio 
of the corrections was taken as unity. 

In the Bethe-Heitler correction term [Equation 2-1] the 
important experimental parameters include the windov; thick- 
nesses and the effective target thicknesses. For this 
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experiment the window thicknesses of the target chambers 
were knov;n (1.0 mil stainless steel for the two-inch dia- 
meter chamber as well as for the three-inch diameter cham- 
ber) . The effective target thicknesses were determined 
from the experimental cross sections, using the assumption 
that the counting system was 100% effective. The calcula- 
tions involved the density of the gas in the target chamber. 
The gas density also appears in the ionization correction 
term. The virial equation of state was used to calculate 
gas densities, to take into account deviations from the ideal 
gas law. The virial equation of state and the method of 
calculating effective target thicknesses are described in 
Appendixes B and C. 

The effective target thicknesses were calculated for 
hydrogen and deuterium, using data from the helium runs and 
Mader's deuterium runs [Ref. 4], The results are shown in 
Table I. 

It was noted from the results that effective target 
thicknesses were somev;hat dependent on energy and very 
dependent on angle. The effective target thickness v;as a 
minimum at 6 = 90°, and increased in proportion to the devi- 
ation of the angle from 90°. For a particular run, the 
effective target thicknesses of the two gases should have 
been equal. Discrepancies of from less than 1% up to 10% 
were noted, however. The deviations v;ere believed to come 
from the margin of error inherent in calculating the areas 
xinder the cross section curves. 
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TABLE I 



EFFECTIVE TARGET THICKNESS 



GAS 


RUN^ 


e 


E]_ (Mev) 


teff (cm) 


H2 


Hel 


90° 


55.96 


.2281.003 


H2 


He2 


90° 


80.10 


.2341.003 




He3 


120° 


85.90 


.8311.012 


H 2 


He4 


120° 


63.18 


.9411.014 


H2 


He 5 


75° 


52.30 


.6831.010 


H2 


001 


90° 


56.82 


.2741.004 


^2 


Okjjl 


0 A 0 


56 . 82 


0 po 4- n n /I 

• ^ <J — • \J \J 'I 


«2 


004 


90° 




.2251.004 


D2 


004 


90° 


85.82 


.2651.004 


H2 


006 


120° 


70.42 


.9001.014 


D2 


006 


120° 




.9441.014 


«2 


009 


120° 


46.55 


.9831.015 


D2 


009 


120° 




1.031.015 



a: RUNS 001, 004, 006 and 009 were H 2 -D 2 mixtures used by 
Mader [Ref. 4]. 
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Gaby [Ref. 3] has used a geometrical approach to calculate 
effective target thicknesses. The geometrical approach, how- 
ever, is very idealized and requires many assumptions, so 
that it is not as reliable as the experimental method. 



B. THE SCHWINGER CORRECTION 



The Schv;inger radiative correction accounts for energy 
loss by electrons in the field of the target nucleus, includ- 
ing nuclear bremsstrahlung , virtual photon exchanges, and 
the emission of photons below the cutoff energy of the 
cross-section calculations. The Schv;inger correction is 
much larger than the Bethe-Heitler correction for the same 
experiment, typically ranging from 10 % to 20 %. 

There are several modifications to the original correc- 
tion developed by Schwinger in 1949 [Ref. 8 ]. The version 
used in NPGLINAC calculations was put forth by Tsai in 1961 
[Ref. 9]. The Tsai version of the Schwinger correction is 




2 a ( 


1/2 In Ej_ + 


1/2 In E 3 


-13' 




In q 2 -1 


- 


n^AE 


AE 


12 




m 2 c 2 



( 2 - 2 ) 




where a = fine structure constant= 1 

137 

= four-momentum transfer squared (in F~^) = q^ - qQ • 
m = electron mass, 

and Ej^,E 3 , AE and n are the same quantities that appeared in 
the Bethe-Heitler correction. The quantity AE was chosen to 
be at least four half-widths of the peak of the cross section. 
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In all the hydrogen-helium runs the half-v/idth was close to 
0.25 MeV, so a value of AE = 1 Mev was used in the calcula- 
tions, not only for the radiative corrections but for 
determining the areas under the cross-section curves as well. 

Tsai also developed a more complicated expression for the 
Schwinger correction.^ It is long and difficult to evaluate 
except by computer. Mader has shown [Ref. 4] that for ener- 
gies obtainable at NPGLINAC, the more complicated Tsai cor- 
rection did not differ significantly from Equation 2-2. In 
fact, for gas mixture e.xperiments where only the ratio of the 
Schwinger corrections for the two gases was needed, the dif- 
ference in the corrections calculated by the complicated 
expression and Equation 2-2 w’as less than 0.1%. Since most 
of the gas target experiments conducted ar NPGLINAC have been 
performed using gas mixtures, it was considered entirely 
adequate to use Equation 2-2 to calculate the Schwinger 
correction terms. 

III. HELIUM FORM FACTOR 

A. INTRODUCTION 

Elastic electron scattering from the He^ nucleus has been 
performed several times, by Hofstadter, et al. [Ref s . 10 , 11] , 
by Erich, Frank, Haas and Prange [Ref. 12], and by Frosch, 
McCarthy, Rand and Yearian [Ref. 6]. At low values of q^ , the 

^ Equation 3-6 in Ref. 5. 
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charge form factor for He^ was found to be in agreement with 
the Gaussian model; 

F (q2) = exp [-<r2>q2/6j (3-1) 

where <r^> is the square of the rms radius of the nuclear 
charge distribution. For He^, <r2>l/2 was determined to be 
1.63 F by Frank, Haas and Prange [Ref. 12]. However, at 
higher values of q2 (q2>6.2 F~2) a deviation from the Gaus- 
sian model was discovered. In fact a diffraction zero was 
found at q2 = 10 F“2 [Ref. 6] . A new expression for the He^ 
form factor v/as obtained on the basis of the experimental 
work of Frosch, et al. 

F(q2) r= jl - (a2q2)i^j exp (-b^/q^) (3-2) 

with n = 6 

a = 0.316 ± 0.001 F 
and b = 0.681 ± 0.002 F. 

The resulting value of <r2>l/2 was 1.68 F. 

For energies obtainable at NPGLINAC, it was only possible 
to study the helium form factor at values of q2 less than 
0.6 F“2. Therefore it was not possible to distinguish betv;een 
the two models. A comparison v/as made to the Gaussian model, 
however, to help determine the efficiency of NPGLINAC with 
respect to gas targets in general. 
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B . THEORY 



Linear accelerators are very useful in producing energe- 
tic electrons which can be scattered from nuclei in order to 
learn something about the charge distributions of the differ- 
ent nucleons and nuclei. The original theories of scattering 
cross sections were developed by Rutherford and Mott, and 
assumed that nuclei were point particles. Nuclei are not 
point particles, however, and their charge distributions have 
a finite size. In the first Born approximation, the ratio of 
an experimental cross section to the Mott cross section 
(including recoil) is called the charge form factor (F) and 
is always less than one for values of q^>0. 

\dQ/ exp/ l.dJ2/ Mott = (3-3) 

where /da \ cos ^2 / 1 \ (3-4) 

\ dJ^ / Mott =\2Eij 0^ \ n / 

sin 4 2 



It can be shown that the charge form factor is just the 
Fourier transform of the charge distribution, and is a 
function of the momentum transfer: 



F(K) = 1 
Be 



jj^xp (-iK *r) p (r) d^r 



(3-5) 



where q is related to K by q = h K 

For a spinless, spherically symmetric nucleus such as He^, 

the charge form factor can be expressed as 

F(r2) =1-1 R2<r2> + ... . (3-6) 

6 
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A plot of F(q ) vs should have a y-intercept of 1.0 and 

a slope of this manner the rms radius of the 

6 

charge distribution can be determined. 

C. EXPERIMENTAL METHOD 
1 . Procedure 

The experiment consisted of using NPGLINAC to scat- 
ter electrons from a hydrogen-helium target at different 
values of q^ and determining the helium charge form factor 
from the resulting data. The values of q^ ranged from 
0.1 F"2 to 0.55 F“2. Appropriate angles and machine energies 
were chosen to give the most reliable data, taking into con- 
sideration the energy limitations of the machine, the length 
of time required for a run, and the background radiation 
levels. 

The data for each run consisted of a spectrum of 
counts per MeV versus electron energy in MeV. Spectrometer 
settings were such that the elastic peaks of hydrogen or 
helium appeared between channels four and five of the ten- 
channel counting system. The total number of points taken 
during a run was dependent on the average energy separation 
between counters. In addition to scanning the helium and 
hydrogen peaks, a number of points were taken above the helium 
peak and between the helium and hydrogen peaks in order to 
obtain measurements of the radiation background levels. The 
integration of the incident beam v/as chosen such that at 
least 10,000 counts were obtained under each peak. The data 
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for each point, including counts for each of the ten chan- 
nels, integration capacitance (pF) and voltage, time of inte- 
gration, incident electron energy (MeV) and spectrometer 
setting (MeV) were recorded on a teletype machine. 

2 . Experimental Errors 

Because the experiment involved calculating the ratios 
of areas under the peaks of helium to hydrogen, systematic 
errors in the experimental parameters cancelled out. These 
errors included the deviation of the incident electron beam 
energy from the set value, which varied less than 0.1% over 
several hours; the drifting of the magnetic field in the 
spectrometer, which was held constant to within 1 part in 
• 10^ by a detailed balancing circuit; and changes in the 
secondary emission monitor (SEM) efficiency. While the ab- 
solute SEM efficiency was not critical for evaluating the 
cross-section ratios, it was more important in determining 
effective target thicknesses of individual gases. 

Statistical errors appeared in both the counting 
rate and the efficiencies of the counters. The efficiencies 
of the counters have been measured by Stewart [Ref. 13]. With 
a counting rate of at least 10,000 counts per peak, the 
statistical error was 1.0% or less. 

The largest error was found to occur when evaluating 
areas under the peaks. This was due to difficulties in 
determining how much of the radiative tail of helium should 
be subtracted from the hydrogen peak. 
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D. DATA REDUCTION 



1 . Calculation of Cross Sections 

The raw data for each run consisted of a teletype 
printout of counts per channel, incident beam energy, spec- 
trometer setting, and integration data. This data was re- 
duced using a computer program which contained counter 
efficiencies and the energies "seen" by each of the counters. 
The reduced data was used to make a plot of energy versus 
counts per microcoulomb. Two examples of these plots are 
shown in Figure 2 and F.igure 3 . They represent values of q- 
of 0.1031 F“^ and 0.5496 F”^ and show the different separa- 
tions of the peaks corresponding to different recoil energies 
of hydrogen and helium. 

The differential cross sections v;ere calculated using 
the following relation: 




(3-7) 



where N 



sc 



number of scattered electrons 



Avg. Energy resolution of counters 



Area under Peak 



X Radiative 
Corrections 



N^ = number of incident electrons 



N^ = number of nuclei per unit area, 
and Afi = solid angle subtended by the spectrometer 



entrance 



It is also known that the experimental cross section equals 
the Mott cross section multiplied by the charge form 
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Scattered Electron Energy (MeV) 

Figure 2. Spectrum of Elastically Scattered Electrons from He - H2 Gas Mixture 
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factor squared lEquation 3-3J . Thus the form factor for a 
gas can be obtained from 



p2 = 



N 



SC 



NiN^AQo^ 



(3-8) 



In this experiment is was desired to calculate the ratio of 
the heliiam charge form factor to the hydrogen form factor. 
Nj^ and AQ fall out of the calculation right away, for they 
are the same for both gases. 



^^'^sc^H ^^t^He ^^M^He 

( 3-9 ) 

The average energy resolution of the counters is directly 
proportional to the energy of the scattering peak. In this 
experiment also, the correction for ionization and the 
Bethe-Heitler correction were the same for both gases. The 
only radiative correction that appeared then was the Schv;inger 
correction Kg. Therefore the ratio of the cross sections 
becomes 



F^He ^ (^^He ^^3^H ^^t^H 

(^3^He (^t^He ^^s^He 

The areas under the peaks were determined using a 
three-point method of integration from a prepared calculator 
program. The appropriate backgrounds were then subtracted 
out. The upper limit on the energy was chosen to be where 
the heliixm curve dropped to background level, which usually 
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5 ? 




was about 0.3 MeV above the He peak. The lower limit was 
chosen to be about four half-widths below the peak energy, 
which was close to 1,0 MeV for all runs. The energy range 
used for the integrations was essentially the same for both 
the hydrogen and helium peaks. 

The gas concentration was calculated from the 
known molar concentrations of H 2 and He. All runs were con- 
ducted using a gas mixture that contained 44.9% hydrogen and 
55.1% helium. There are two H 2 nuclei for each He nucleus. 
Thus the ratio of nuclei concentrations was 

(N^.)„ ^ 2(44.9) ^ 1.62976 

(Nt>He (55.1) 

The Mott cross sections were computed from 
Equation 3-4 using the desk calculator, and the Schwinger 
radiative corrections were calculated from Equation 2-2. The 
hydrogen form factors were obtained from a prepared 
calculator program. Finally, the helium charge form factors 
were calculated from Equation 3-10. 

2 . Error Analysis 

All errors in the calculation of the helium charge 
form factor were 1.0% or less except for the error in eval- 
uating areas under the cross-section curves. Specifically, 
it was difficult to determine exactly how much to subtract 
out from the hydrogen peak due to background and the helium 
radiative tail. One method used v;as to take an average back- 
ground rate between the two peaks and subtract it frcm the 
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hydrogen peak. At low values of , however, the peaks were 
very close together and this made it difficult to determine 
an actual background rate [Figure 2] . Another method used 
was to make a least-squares fit of the helium tail from 
halfway up the peak to the start of the hydrogen peak. A 
program was prepared to perform a least-squares fit of the 
form f(x) = a + bx + cx^ , 

where Xj_ = 1 , Eq being the helium peak energy, 

Eq - 

and f(xi) = counts per microcoulomb corresponding to energy 
E^. The least-squares fit v/orked well at low values of q^ , 
but at q2 = 0.55 F“^ the quadratic terra caused the fit to 
begin rising under the hydrogen peak. In this case the fit 
was reduced to f(x) = a + bx and results were better. 

The differences in calculations using an average 
background rate and a least-squares fit were up to 10%. The 
method of least squares resulted in form factor values more 
closely correlated to the expected results. Therefore the 
least-squares method was used for all calculations in this 
experiment. 

E. RESULTS AND CONCLUSIONS 

2 2 

The experimental values of F(q ) for each value of q 
used are listed in Table II and are depicted graphically in 
Figure 4. The error flags indicate the maximum statistical 
error of 1.5%. The theoretical curve for a Gaussian poten- 
tial is also shown for comparison. Although the experimental 
values of F(q ) differed somev/hat from the theoretical values 
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TABLE II 



EXPERIMENTAL He^ CHARGE FORM FACTOR 



RUN 


q2 (f~2) 


(Fh^)- exp 


(Fhq) Gaussian 


1 


.1585 


. 982± . 015 


.9322 


2 


. 3328 


. 813± . 012 


. 8630 


3 


. 5496 


. 783±. 012 


. 7840 


4 


. 3000 


. 815±. 012 


. 8756 


5 


.1031 . 


. 957±. 014 


.9554 
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Figure 4. He^ Charge Form Factor vs. 



at low values of , a least-squares fit to the form 
y = ae^^ produced a curve that nearly coincided with the 
theoretical curve. The resulting rms radius of the He'^ 
nucleus, obtained from the least-squares fit, was <r2>l/2 
= 1.78 F, larger than the expected value by 9.2%. The dis- 
crepancies were believed to be due mainly to the difficul- 
ties encountered when evaluating areas under the hydrogen 
peaks. At low values of q2 the separation between peaks v;as 
not wide enough to obtain an accurate determination of back- 
ground and helium radiative tails. It was noted also that 
for q2 = .1585 F“2 the background level increased signifi- 
cantly from the helium peak to the hydrogen peak, causing the 
experimental form factor to be larger than predicted. At 
q = .53 however, cxiere was good agreement between experiment 
and. theory. This vr/as believed due to adequate separation 
between peaks for an accurate background determination, and 
to a fairly constant background rate throughout the energy 
range of the calculations. 

It is believed that better agreement between helium 
form-factor experiments at NPGLINAC and those at other lab- 
oratories can be obtained by using pure gas targets instead 
of mixtures. This v7ould ease the burden of attempting to 
calculate the proper background rates. 



29 



APPENDIX A 



FREQUENTLY USED CONSTANTS 

The following physical constants appeared several 
in cross-section calculations and are listed here for 
reference : 

e^ = 0.144 X 10”^^ MeV - cm 
(Mc2)jj = 938.68 MeV 
(Mc2)He = 3728.01 MeV 
h = 

ra^c 386.12 F 

2a = _ 

~ 0.46455 X 10“^ 

] cm^ = 

erg 0.68947 x 10^ in.^ 

lb 

RT = 0.643041 x 10^° erg at T = 77.32° K 

mole 

1 d yne = 0.14504 x 10"^ lb 

in . 2 



times 

future 
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APPENDIX B 



VIRIAL COEFFICIENTS 

The theoretical development of the virial theorem and 
the virial equation of state is described in Ref. 14 and 
Ref. 15. The virial equation of state is an alteration to 
the ideal gas law (pV = nRT) and more accurately describes 
the properties of real gases. The virial equation of state 
appears in two forms, a power series expansion in the molar 
volume and a power series expansion in the pressure. Since 
pressure is a directly measureable quantity for gas target 
runs at NPGLINAC, the latter form is more convenient to use: 

£V - 1 + (T)p + C" (T)p2 + . , . . (B-1) 

RT 

V = V , the molar volume, where Nq is Avogadro ' s numbeir. 

No 

B"(T) = B/RT and C" (T) = (C-B^) / (RT)^, 

where B and C are the second and third virial coefficients 
and are functions of temperature only. For the gases current- 
ly being studied at Monterey (hydrogen, deuterium and helium) 
the virial coefficients were calculated from the Lennard-Jones 
potential for non-polar molecules. Preliminary calculations 
showed that for a pure gas target, the values of B'’(T) were 
of the order of 10 in. /lb and values of C^(T) were of the 
order of 10 in. /lb. The terms involving higher orders of 
p were dropped with no significant changes in results. 
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Evaluation of B(T) and C(T) for pure gases is straight- 
forward, using the relations 

B = B*(T*) and C = C* (T*) . 

^ 0 

bQ = 2/3 ttN|^ cm^/mole, and T* = kT/e, 

the reduced temperature, a, e/k and bQ are force constants 
for the Lennard- Jones potential and are tabulated for sever- 
al gases [Ref. 15], including those in use at Monterey. Two 
values for each parameter are quoted, one determined classi- 
cally and the other quantum mechanically. The quantum- 
mechanical parameters are more accurate in principle and were 
used in all calculations. 

Once T* has been determined, B* and C* can be extracted 

from a table [Ref. 15]. For values of T* not exactly equal 

to the tabulated values of B* and C* , a linear interpolation 

is necessary to obtain correct values of B* and C*. Proper 

substitution yields the quantities B^(T) and (T) , v.’hich can 

then be entered into the equation of state. All parameters 

are then known except V, v;hich is related to the desired 

parameter p (density) by V = W, where W is the molecular 

P 

weight. Solving for the density yields the following 
formula : 



P 



I = ^ 

^ RT [1 + B"(T)p + C"(T)p2] 



(B-2) 
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The calculations of the virial coefficients and 



Cj^ix for a gas mixture are slightly more complicated. In fact, 
C^ix can not be calculated exactly so it is dropped from the 
expansion. In general, f°^ ^ mixture is given by 

n n 



B 



mix 




Xi Xj Bij , 



(B-3) 



i=l j=l 



where n = number of constituent gases in the mixture and 
= molar fraction of the ith gas in the mixture. 

B^j (i ^ j) is a theoretical virial coefficient for the gas 
mixture, v;here Tij* arid bQ are determined from the new para- 
meters = /FT FT and = 1/2 (a + Oj) , obtained from 

the force constants of the ith and jth gases in the mixture. 
Bj^j* is extracted from the table in the same manner as be- 
fore, only the table is entered with the value 
Tij* = kT/Sj^j . Bjj is merely the second virial coefficient 
for the pure jth gas in the mixture. 

For typically used gas mixtures (n = 2) , the second 
virial coefficient becomes 

Bmix = + 2Xj, Xj (B-4) 

The temperature of all gases and gas mixtures used at 
NPGLINAC was the temperature of the liquid nitrogen that was 
used as a coolant (77.32° K) . Variations from this tempera- 
ture were small during all runs, and it was assumed constant 
for all calculations. Figure 5 depicts the temperature 
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4 






bridge control mechanism. Values of the virial coefficients 
for the gases used are listed in Table III. Table IV shov;s 
values of gas densities for typical pressures calculated 
from the virial equation of state. 
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TABLE III 



VIRIAL COEFFICIENTS FOR SEVERAL GASES AT T = 77.32° K 



j. 



^ 

GAS 


B"(T) (in.2/ib) 


C" (T) (in.'^/lb^) 




-.18966 X 10"^ 


+ .12928 X 10"’^ 


^2 


-.18966 X 10 “^ 


+ .12928 X lO""^ 


48.8% H2~ 
51.2% D 2 
mixture 


-.18944 X 10"^ 




\ 

He 


+ . 89735 X 10“"^ 


+.71852 X 10“^ 


44.9% H2“ 
55.1% He 
mixture 


+.43683 X 10"5 
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TABLE IV 



DENSITIES OF SEVERAL GASES AT CONSTANT 
TEMPERATURE (77.32° K) FOR VARIOUS PRESSURES 



Gas 


p(lb/in. 2 ) 


p(gm/cm^) x 10 ^ 


H 2 (in 


140 


. 310906 




148 


. 329184 




149 


. 331473 


H 2 -D 9 


150 


. 333699 




151 


. 336053 


mixture) 


152 


. 338344 




160 


. 356709 




140 


. 621334 


D 2 (in 


148 


. 657862 




149 


. 662437 


H 2 -D 2 


150 


. 666885 




151 


.671591 


mixture) 


152 


.676170 




160 


. 712871 




140 


. 302475 


H 2 (in 


148 


. 319748 




149 


. 321907 


H 2 ~He 


150 


. 324066 




151 


. 326225 


mixture) 


152 


. 328384 




160 


. 345655 




140 


. 600643 


He ( in 


148 


.634943 




149 


. 639231 


H,-He 


150 


. 643518 




151 


. 647805 


mixture) 


152 


.652093 




160 


. 686389 
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TABLE IV (Cont.) 



DENSITIES OF SEVERAL GASES AT CONSTANT 
TEMPERATURE (77.32° K) FOR VARIOUS PRESSURES 



Gss 


p (Ib/in. ^) 


p(gm/cm^) x 10 ^ 




140 


. 310834 




148 


. 329099 




149. 


.331386 


(pure) 


150 


. 333839 




151 


. 335963 




152 


. 338252 




160 


. 356600 




140 


. 621204 


^2 


148 


. 657694 




. 149 


.662264 


(pure) 


150 


. 667166 




151 


. 671411 




152 


. 675986 




160 


. 712654 




140 


. 593471 


He 


14 8 


. 626929 




149 


. 631108 


(pure) 


150 


.635286 




151 


. 639463 




152 


. 643640 




160 


. 677024 
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APPENDIX C 



GAS TARGET EFFECTIVE THICKNESS 



The thickness of a solid target is directly measureable 
and poses no problem in a cross-section calculation. How- 
ever, a gas target is different. The pressurized gas is con- 
tained in a cylindrical chamber two or three inches in . 
diameter, but only those electrons scattered from a small 
volume of the gas make their way through the spectrometer 
entrance slit. The calculation of the effective thickness 
of a gas target is detailed, and is reviewed here to show 
what was done for this experiment, and to clarify the 
procedure for future reference. 

The experimental cross section is given by Equation 3-7. 
The target thickness t is located in the term N^, the number 
of target nuclei per unit area. 

Nt = ptNQ (C-1) 



where A is the atomic weight, p is the gas density determined 
using the virial equation of state, and Nq is Avogadro's 
number. Substituting Equation C-1 into Equation 3-7 and 
solving for t gives the following expression: 



A 



t 



N 



sc 



A 




(C-2) 
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I 



5 % 



t 




Ng^ is the number of scattered electrons and is equal to 
the integrated area (counts-MeV) divided by the average reso- 
lution of the counters AE^. (MeV). . AE^ can be calculated from 
the dispersion formula, 

AE^ = 

E p D (C-3) 

where a = 3/16" 
r = 16" 

and D = Dispersion factor = 3.92. 

Thus E^= 0.00299 E. 

The solid angle AQ = Area of spectrometer entrance slit/r^ 

= (5/16) (1 1/2) 

162 

= 0.183 X 10“2 

The number of incident electrons is determined from the total 
charge accumulated in the SEM and the SEM efficiency: 

Ni = Q 

e (C-4) 

e is the charge of one electron. The SEM efficiency was 
assumed constant at 6.1%. 

Now the actual integrations calculated from the reduced 
data were in units of counts per yC - MeV. This causes the 
Q in Nj^ to drop out when N^^/Nj^ is calculated: 

Ngc = Area-e-EgEM 

Nj_ AEj. (C-5) 
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The experimental cross section is evaluated using Equation 
3-3. Thus the final form of the expression for effective 
target thickness becomes: 

t = Area*e*E5£j.|*A 
*AE^*Ar2*p *Nq 

= (0.29605 X 10 ~ 32 ) f Area • a \ cm. (C-6) 

\ / 
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